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We calculate total hadron-hadron cross sections at high energies using the Low-Nussinov two- 
gluon model of the Pomeron. The gluon exchange is represented by a phenomenological potential 
including screened color-Coulomb, screened confining, and spin-spin interactions. We use bound- 
state wave functions obtained from a potential model for mesons, and a Gaussian wave function for 
a proton. We evaluate total cross sections for collisions involving n, K, p, (f>, D, J/ip, ip', T, T', 
and the proton. We find that the total cross sections increase with the square of the sum of the 

C*") ' root-mean square radii of the colliding hadrons, but there are variations arising from the spin-spin 

interaction. We also calculate the total cross sections of a mixed-color charmonium state on a pion 
■ and on a proton. The dependence of the total cross section on the size and the color state of the 

s ! ' charmonium is investigated. 
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PACS numbers: 13.85.Lg, 14.40.-n, 14.20.-c, 12.39.- 



I. INTRODUCTION 



Hadron-hadron total cross sections are important characteristics of the dynamics of hadronic systems. Theoretical 
descriptions of the reaction process, however, remain incomplete as the mechanism of how hadrons interact contains 
important aspects of non-perturbative QCD. A phenomenological description of Donnachie and Landshoff Q] gives the 
total cross section in terms of the exchange of a Pomeron and a Reggeon. The exchange of a Pomeron dominates the 
cross section at high energies. It corresponds to the exchange of a tower of fictitious particles on a Regge trajectory 
with a Regge intercept of approximately 1. The exchange of a Reggeon dominates the cross section at low energies. 
It corresponds to the exchange of p, u>, fi, a^... on the Regge trajectory with a Regge intercept of approximately 0.5. 

It is desirable to study microscopic models of total hadron-hadron cross sections at high energies. The total cross 
sections provide useful information on the sizes of hadrons and the dynamics between them. Furthermore, they 
are important input data for the investigation of other reaction processes. For example, in high-energy heavy-ion 
collisions, the suppression of J/ip production has been proposed as a signature for the production of the quark-gluon 
£h ■ plasma 2]. However, the suppression of J/tp production can also arise from the collision of J/ip with hadrons at high 
energies. It is necessary to understand the absorption of Jhp by hadrons before one can unambi guo usly identify the 
quark-gluon plasma as the source of J/tp suppression [lllHll0HI3[l3[ll0mEl ^ ^ El El El HJ • 
While the systematic studies of Donnachie and Landshoff provide very valuable information on the total cross section 
of many hadron-hadron systems, there are reactions (involving, for example, J/ip) for which data are not available. 
They can only be evaluated theoretically |2l], \'2'j l23t l24l 123. 12(1 127| . It is therefore useful to develop a microscopic 
theory of hadron-hadron collisions at high energies. 

As the Pomeron exchange is the dominant process in high-energy hadron-hadron reactions, it is useful to model the 
exchange of a Pomeron explicitly in microscopic terms. Previously, the Pomeron exchange in a hadron-hadron reaction 
was studied in terms of the exchange of two gluons. The two-gluon model of the Pomeron (TGMP), first proposed 
by Low [2^| and Nussinov [2^, was further investigated by Gunion and Soper [3(j who introduced an effective gluon 
mass to mimic the confinement of the colored gluon. They examined the effects of the size of the bound-state hadrons 
and the number of quarks they contain. Landshoff and Nachtmann generalized the concept of the gluon condensate 
to the correlation length of the gluon field in the vacuum and discussed the short-range nature of the effective gluon 
exchange between quarks in hadron-hadron scattering |3lJ . The dependence of the cross sections on the sizes and the 
color states of the colliding hadrons was considered by Dolejsi and Hiifner [3i| • They obtained useful analytical and 
numerical results for the total cross sections of many hadron-hadron systems. The effects of channel coupling for the 
propagation of a radially-excited hadron |3^ | was investigated in the two-gluon model of the Pomeron by Wong [34| 
who also evaluated many meson-nucleon total cross sections. 

We shall follow the approach of Gunion, Soper, Landshoff, Nachtmann, Dolejsi, Hiifner, and Wong to study the 
hadron-hadron total cross section, with the following refinements. We shall use a set of meson wave functions obtained 
from meson mass calculations [lg to replace the simple form factors used in earlier studies. This allows a more 
detailed and systematic study of the hadron-hadron scattering processes involving a much large set of hadrons than 
those investigated previously. Gunion and Soper |3(j, Dolejsi and Hiifner :32], and Wong [34[ made use of only a single 
screened color-Coulomb potential to represent the non-perturbative gluon exchange. We shall instead represent the 
gluon exchange between the constituents of one hadron and constituents of the other hadron by a phenomenological 
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potential containing screened color-Coulomb, screened color-confining, and spin-spin interactions |l8l 13a . 1361 l37| . 
Screening arises when the interaction between a constituent of one hadron and a constituent of another hadron 
occurs at large distances for which the production of dynamical quarks screens the interaction. The use of a more 
complicated potential allows a greater degree of flexibility, which provides a better description of non-perturbative 
effects in scattering. It also enables one to use reasonable strong interaction coupling constants, as the constraint 
of using a single screened color-Coulomb interaction may sometimes lead to values of a s larger than unity [32j. In 
addition, the spin-spin interaction is known to be important in non-relativistic quark models. Its phenomenological 
incorporation here will allow us to evaluate the variation of the total cross section in reactions of mesons with different 
total spin quantum numbers and quarks of different flavors. 

Similar to earlier works of Dolejsi and Hiifner [3^, and Wong [34j], we shall first obtain the the screening mass and 
the strength of the interaction to reproduce the experimental pp, np, and Kp cross sections at y/s — 20 GeV. We then 
use the theoretical meson wave functions to calculate other cross sections. Based on the energy dependence of the 
phenomenological description of total hadron-hadron cross sections found by Donnachie and Landshoff 0] , we shall 
then extend our results to higher energies. 

In reactions involving heavy quarkonium production in nuclear collisions, the initial heavy quark pair produced as 
a result of a parton-parton collision is a coherent state of a color admixture and different projections of this state 
will lead to different final heavy quarkonia 0, |3^, 13^. liol Ell l42l |43| . We will also examine the color dependence 
of the cross section, in order to understand how a colored heavy quark pair propagates in a hadronic medium 
[H El3- EM DH EH HH E3 • The nucleons with which the produced charmonium collide can also acquire a color as a 
result of prior collisions and may collide with a color-singlet charmonium [42l l43j| . In this paper, we shall limit our 
discussion to the problem of the interaction of a colored charmonium state with a color-singlet hadron. 

Our approach to study the total cross section in the TGMP differs from that of the additive quark model 
and the dipole model of photoproduction ^5j. The total cross sections obtained in different models can naturally 
be different. Careful comparison of these different results with experimental data will be useful to determine the 
importance of various scattering mechanisms that are present in the interaction of hadrons. 

This paper is organized as follows. In Sec. II we describe the model used to calculate the hadron-hadron elastic 
scattering amplitude at high energies. In Sec. Ill we evaluate the spin matrix element for meson-meson, meson-proton, 
and proton-proton scattering. The calculation of the spatial matrix element is discussed in Sec. IV. In Sec. V we 
present results for total hadron-hadron cross sections and the energy dependence of the total hadron-hadron cross 
sections. In Sec. VI we evaluate the total cross sections of a mixed-color charmonium state scattering on a pion or 
a proton. The dependences of the total cross sections on the size and the color of the initial charmonium state are 
investigated. Finally, we present our summary and conclusions in Sec. VII. 



II. TWO-GLUON MODEL OF THE POMERON 



To calculate the total cross section in a hadron-hadron collision, we consider the hadron-hadron elastic scattering 
amplitude A. According to the optical theorem, the imaginary part of the elastic scattering amplitude at zero 
momentum transfer gives the total hadron-hadron cross section. 

The elastic scattering process between two hadrons at high energies is dominated by the exchange of a Pomeron, 
as evidenced by the slow variation of the total cross section with energy, the property of having an approximately 
purely imaginary forward scattering amplitude, and the dominance of no quantum number flow in the fragmentation 
region in forward inelastic processes. The phenomenological treatment of the total cross section as arising from 
the exchange of a Pomeron with a Regge trajectory intercept of 1.0808, and the exchange of a Reggeon from the 
trajectory of PjJjJ, fi, a^... with an intercept of 0.5475, gives an excellent representation of many hadron-hadron total 
cross sections [l,]. The small deviations of these Regge trajectory intercepts from the lowest-order expectations of 1.0 
and 0.5 respectively have been attributed to higher-order effects 0. 

In this work, we shall model the exchange of a Pomeron in terms of the exchange of gluons. As the exchange of 
a single gluon would lead to an exchange of color, the lowest number of gluon exchanges with no net color flow is 
the exchange of two gluons, proposed first by Low |2§] and Nussinov [2£j and studied by Gunion, Soper, Landshoff, 
Nachtmann, Dolejsi, Hiifner, and Wong |30L l3lL 133. l34f . 

Accordingly, we express the total hadron-hadron cross section in terms of the hadron-hadron elastic scattering 
amplitude at zero four-momentum transfer squared t 

cr tot = -Im A(s,t = 0). (1) 
s 

Following Gunion and Soper [30l | and Dolejsi and Hiifner |32| . we obtain the forward scattering amplitude A by 
including diagrams of the type shown in Fig. 1. As the elastic scattering process does not change the hadron internal 
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structure, the initial and final hadron states are identical. In Fig. 1, constituents i or j in hadron I interact with 
constituents k and I in hadron II by exchanging two gluons. As the elastic process involves important aspects of 
non-perturbative QCD, it is reasonable to represent the corresponding gluon exchange between constituents i and k 
in terms of a phenomenological potential V(ik) possessing non-perturbative elements of color-Coulomb, confinement, 
and spin-spin interactions. The complete amplitude is obtained by summing over all possible diagrams representing 
the exchange of two gluons between all constituents {i, j} of hadron I with all constituents {k, 1} of hadron II [32|: 



A(s,t) 




d 2 b 



,-iQb 



i,j,k,l 



\CiC llXl Xn *I *n| V T (i*0 



y-V T {jl)\CiC u Xl Xi *i *u! 



(2) 



where b is the impact parameter, Q 2 = — £, C\ and C\\ are the color wave functions, X \ an d X u are the spin wave 
functions, and and \I>n are the spatial wave functions of hadron I and hadron II, respectively. The square brackets 
represent the averaging over the color, spin, and internal spatial degrees of freedom of hadrons I and II. The sum 
J2i j k i ^ s over au two-gluon exchange diagrams. 

At high energies when the average intrinsic quark momentum inside a hadron is small compared to the center- 
of-mass energy, the case of a small momentum transfer in a two-gluon exchange diagram can be approximated as a 
loop diagram in which the quark lines on the loop can be treated as on the mass shell and represented by two delta 
functions (as in Eq. (A. 11) of Ref. HJ). After integrating over the time-like and the longitudinal loop momenta using 
these two delta functions, the resultant effective gluon propagators involve only transverse momenta (or equivalcntly 
only transverse coordinates) [30L l3lL 132^ . When we represent the effective gluon propagator phenomenologically by 
an interaction, the interaction Vr(ik) in transverse coordinates, written explicitly as Vt(i"t), can be obtained by 
integrating the interaction V(r) over the longitudinal coordinate z, 



V T {ik) = Vr(r T ) 



dz V(r) 



(3) 



where r is the three-dimensional relative coordinate between the interacting particles i and k and rx its transverse 
component. Defining the Fourier transform of V(r) as t^(k) by 



V(r) = 



dk 



we can obtain from Eqs. © and Q 



V T (r T ) 



dk 



T „iki 



(2tt) 2 



V(k) 



(4) 



(5) 



k z =0 



Thus, Vr(kx), the 2-dimensional Fourier transform of Vt(i"t), is just the 3-dimensional Fourier transform V"(k) 
evaluated at k z — 0. 

As the interaction between the constituents takes place at large distances, we are well advised to use a screened 
potential to represent the effects of dynamical quarks. As pointed out by Landshoff and Nachtmann |3l| . two quarks in 
two hadrons can exchange a nonperturbative gluon only if they pass within a short distance of each other. This short 
correlation length of the gluon field is another manifestation of the screening phenomenon . We employ screened color- 
Coulomb (Yukawa) and screened confining (exponential) potentials [35l 13a . l37| and a spin-spin interaction 0, [3^ 
for our interaction V(r). The interaction between quark i in hadron I and quark k in hadron II in three-dimensional 
relative coordinate r can be written as 



V(r) 
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where is the Gell-Mann matrix, and is replaced by — X i for an antiquark i. The quantities a s , \i, bo, and 
d are the strong coupling constant, the effective screening parameter, the effective string-tension, and the spin-spin 
potential width parameter. The Fourier transform of the spatial and spin part of the potential is 



V(k) = 



k 2 +/i 2 (k 2 + i u 2 ) 2 



87ra s 
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3mimfe 



s k e 



-k 2 /4d 2 



(7) 
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We have used an exponential interaction in Eq. JfJJl to represent the screened confining interaction as it contains the 
appropriate properties. At short distances for which screening effect is small, the exponential interaction gives the 
linear interaction b$r, and at large distances for which the confining interaction is screened by the production of 
dynamical quarks, it approaches zero. In terms of the Fourier transform, screening is represented by a non-vanishing 
value of the screening mass [i is Eq. JjJ. The use of a screened confining interaction is necessary here. If we assume 
an unscreened linear-confining interaction, corresponding to taking /j = in the 6-7r&o/(kT + /i 2 ) 2 term in Eq. Q), 
we shall see in Section V that the total hadron-hadron cross section will be singular. 
The square brackets in Eq. @ can be written as 

; Cl C„ X, Xu *I *H I M^) V T (jl) I d C U Xl X: *I *II ] 
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(c ijkl J2(Ci\ XfX'j I Ci ){ Cii | X a k \\ \ C u )}{ ^ E 



x (x : Xu *i(ri) *n(r n ) I V T (r lkT )V T (r jlT ) \ Xl X„ *i(ri) *n(r n ) )}, (8) 



where the first pair of curly brackets is the sum of color matrix elements, Cijki is the product of the four factors Ci, 
Cj, Ck, and Ci, each of which is 1 for a quark and —1 for an antiquark, and a and b are gluon color component labels. 
In the second pair of curly brackets, the sum is over all spin states of hadrons I and II (total number N x ). 
The color wave function for a color-singlet meson is 



C 1 ) = -^^|aa), (9) 

* a—l 



where a (a = 1,2,3) is the quark color label and a is the corresponding antiquark color label. The color matrix 
element is 

(C 1 |A^| C 1 )-^ 6 . (10) 

For a color-singlet nucleon we have 

1 3 

l Cl ) = ^ E e a/97 |a/3 7 ), (11) 



V6 

a,/3,7=l 



and 




|a;„ a!;i r" ={'•; ,;, f r m " n ' (12) 

tor wiytfi. 

III. SPIN MATRIX ELEMENTS 

From Eqs. and ©, the matrix element for the spatial and spin part can be separated into a sum of products 
of spin matrix elements and spatial matrix elements. The first-order spin-spin matrix elements and the second-order 
spin-spin matrix elements are 

4k = {Xi X„ I (Si • s fc ) I XiXn), (13) 



4k]i = ( Xi Xu I (Sj • B k )(sj ■ si) I Xi Xu > • (14) 
We shall give results for these quantities for meson-meson, meson-proton, and proton-proton collisions separately. 
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A. The meson meson case 

Consider meson I with constituents i and j coupled to total spin Si, and meson II with constituents k and I coupled 
to total spin Su. We would like to recouple the constituents so that the operators (s, • Sfe) and (sj ■ s;) have simple 
eigenvalues. We have 

I Xi Xn ) = I (klfsll ) = 52(SS z \SiSi z SnSn z ) | [(ij) s >(kl) s »}l > 

= E E (SS^SiSi^SiiSii^SiSuSikSjt I si s\ Su 1 | {(ik) s **(jl) s »}l ) 

= E E (SS^SiS^SiiSii^i-lf^-^SiSiiSuSjkl s\ Z Su 
s,SzS u ,s jk [Su Sjk S 

x\[(il) s Hjk) s **f Sz ), (15) 
where S = ^25+1 . From the above equation, the matrix element s^} — (\i Xu I ( s « ' s *=) I Xi Xu ) IS 

4k =E E liSSzlSiSizSuSii^SfSlStJ^lZ t Iil I[5 ifc (S ifc + l)-|], (16) 



s,SzS ik ,Sji [S ik Sjt S 

and the matrix clement sf^ = ( x, Xn I (s» ■ s fe )(sj ■ s/) | x, Xu ) are 



s (2) 

ikik 



E E i(^ift^z^ii z )i 2 ^i^^| S 5 til {\[Sik{s ik + i)-h\ , (l?) 

S,S,S,>,S„ [S ife S jt S J ^ j 

E E E i(^^i^iz^ii5„ z )i 2 (-i) Sii - sfc - s '(-i)^ ! - s ^ s '(-i)^ fc - s ^ sfc 

( Si s-j Si] ( Si s k Sik} [ Si sj Si 
xSjS n SlS 2 JlS%{ s k Sl SnSlst S j Sjt} 1st s k Su 



Sik Sji S J I Su Sjk S J I Su Sjk s 

2 J " 2 



x \ [S tk {S lk + !)-§] x \ [Su(Su + !)-§], for fc ^ / , (18) 



s,s z s ik ,s } i [S ik Sjt S J 

x^[^(5,/ + l)-|, fori^j andfc^Z , (19) 



where i, j = 1,2; k,l = 3,4. 



B. The meson proton and proton proton cases 



The spin and flavor wave function for a proton in the S z — 1/2 state is 

|Pi/2> = y= ( 2|u+u+eL) - \u+u-d+) - \u-u+d+] 



+ 2 | u + d-U + ) — | u + d + U- ) — | U-d + u + ) 

+ 2\d-U + u + ) - | d + u + u- ) - | d + d-u + ) ) . (20) 
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Denoting the two it-quarks by 3 and 4, the d-quark by 5, and the spinor of quark i as i s s , we have 



1/2' 

h/2 • 



( 2 I 3 



!/2 ,1/2 (.1/2 , 
1/2 *l/2 -1/2 / 

1/2 rV2 .1/2. 



Z I °l/2 °-l/2 *l/2 

2 I 5 1 / 2 3I/2 4 i/2 v _ I g i/a 3 i/2 4 i/a 



ql/2 A l/2 t-1/2. 
J l/2 *-l/2 °l/2 / 

ql/2 k1/2 ,1/2 , 
J l/2 °l/2 *-l/2 1 



ql/2 ,1/2 ,1/2, 
°-l/2 *l/2 °l/2 1 

ql/2 cl/2.1/2, 
°-l/2 d l/2 n/2 ' 



I cl/2 ql/2 4 l/2n 

•"-1/2 "1/2 ^1/2 / I "1/2 "1/2 ~*-l/2 / I °l/2 -1/2 1/2 / ) ' 

We can couple the spin of any two of 3, 4, and 5. For example, coupling 3 and 4, we have 
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1/2 . 
1/2. 
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E 
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2 2 2 2 
/111 
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V 2 2 2 



(34)1 5^ /2 ) (1) + | (34)} 5%) w + | (34)} 5^ /2 ) (7) 

S34 rl/2 \ l/'Q/l\534 rl/2 



34- 



l(34)^5^) (2) - |(34)g- 5^) w 



|(34) 



S34 rl/2\ 

S 3 4, °l/2/(8) 



E 

S34,5 , 34 2 



|(34)5) 



1 1 1 

222 



S34S34Z 



\(M)kA/l)&)- 1(34) 



S34 cl/2 
S 3 4z °l/2 



)(6) 



|(34) 



1st 5 l/2)(9) 



(21) 



(22) 



where the subscript (n) labels the flavor component in the n-th term in Eq. I|21|) (i.e., („)( (W)s*J ' I (^)s'' ' /(«•') = 

8{SkizS' klz )8(SkiS' kl )5{nn r ) ). We can also write out the other coupling expressions of | p\i^), for example, | (35)4) 

and I (45)3). One can obtain a similar decomposition for | p_jy 2 )• 

Thus, the first-order and second-order spin-spin matrix elements for a meson-proton scattering can be expressed as 

4l ] = <fo-)!i m)i']\(si-s k )\(ij) s s i [(H) 1']) 

= ((ii)f,(w)l(*i-«fc)l(y)S(«)>. (23) 



s (2) _ 



<(W [(fcOni(*<-*iO(*i-*i)l(OT, P0H> 
( («) I (*i-«*)(«i-*0 1 («)£(«)). 



(24) 



where /' is a spectator, j,j = 1,2, and k,l = 3,4,5. Each term in Eqs. i12.'it and 124J1 can be calculated in the same 
way as in the meson-meson case. For a proton-proton scattering, we have 



,(i) _ 



<mf] w)i'}\(si-s k )\[mj'} w)n 



= ((ij)(kl)\( Si -s k )\(ij)(kl)), 



(25) 



= ((ii)(W)l(*i-*fc)(*i-*i)l(y)(«)>. 



(26) 



where f and /' are spectators, i,j = 1,2,3, and k, I — 4,5,6. Each term in Eqs. i|25J) and (|26J) can be calculated 
similar to the meson-meson case. 



IV. SPATIAL MATRIX ELEMENTS 



A. The meson meson case 



Using a set of Gaussian basis states with different widths, Wong et al. previously obtained bound-state meson wave 
functions from a non-relativistic quark model |l8j . The model assumes color- Coulomb, linear-confining, and spin-spin 
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interactions. It gives meson masses in reasonable agreement with experimental data. The meson wave functions are 
tabulated in Table IV of Ref. 0]. We shall use these meson wave functions in our calculations. They are represented 
by a nonorthogonal sum of Gaussian basis functions of different widths, 



,n/3 2 , 3/4 _i3 

n— 1 n—1 



*(r) = J>*,(r) = £ a n (^-) 3/ V^ 2 , (27) 



where r = (rx, 2) is the three-dimensional coordinate, a n is the coefficient of the n-th component, and j3 is a parameter 
characterizing the width of the basis functions in momentum space. 

The transverse coordinates of the quarks and antiquaries in meson I and II can be expressed as 

b riT b r nT roQ , 

riT,2T = — ± — , r 3Ti 4 T = -— ± — — , (28J 

where rix or tut is the relative transverse coordinate of the quark and antiquark in meson I or II respectively. The 
relative transverse coordinate of quark (or antiquark) i in meson I and quark (or antiquark) k in meson II can be 
expressed as 

rifcT = r iT - r kT = h + — n T + — ~- niT , 

(* = 1,2, fc = 3,4, C n = l, C I2 = -1, C II3 = -1, C„ 4 = l.) (29) 
From Eqs. (0), JSJ, (JSJ), l|27(l . and (|29(l . the spatial matrix element for a meson- meson scattering is 

J d 2 b e- jQ b /" dn dr„ *J, f (ri)*J [if (rn)V«r(r ftT )V«r(r J -rr)*i4(ri)*n,i(rn) 

= 7^2 / d 2 k T V % T (k T )VT(k T -Q)^(Q,k T )F fc 7(Q,k T ), (t,i = l,2;*,l = 3,4), (30) 
where the superscript in F^f(Q) labels the colliding hadron as a meson M or proton P, and 

^(Q,k T ) = £a„ a „^^j ex P | 8(n + n ,)ff j. ( 31 ) 



ri . n 



pM , n , , V- /2VW\ 3/2 f [k T + g IIfc C m (Q-k T )] 2 I 

^fcz (Qi k T = 2^ a m a m > exp <^ f . 32 



m,m' 



B. The meson— proton and proton— proton case 

The three-dimensional coordinate of a quark in a proton can be written in terms of Jacobi coordinates r\ = 
(VttVz) an d £ = (Cti&s) P^ - Its transverse component is given by 

r feT = -^ + -r»7 T -|-— -| T , (33) 



2 6 11 2 



where 

[2, fc = 3, 



-1, fc = 4, 5, 



fc = 3, 

Ou : < I. k = 4, (34) 

jfc = 5. 



The relative transverse coordinate of the quark (or antiquark) i in meson I and the quark k in proton II is 

rikT = r iT - r k T = b + — rix — T7 X — £ T . (35) 
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We use a simple Gaussian wave function for a proton. In terms of Jacobi coordinates, it can be expressed as [32 



*(tj, = ( 24V3 r 2 ) 3 T 1 / 2 exp[(-7, 2 /12 - £ 2 /4)/( r 2 p ) ] . 
From Eqs. J3J), ©, JSJl, (|27|l . and l|34f> — l|36[l . the spatial matrix element for a meson-proton scattering is 

J d 2 b e- iQ b / dri dry d£ tf^ (n) *l 1<{ (ri,£,)V T {r ikT )V T {r jlT ) ^(n^ifa, £) 

= (2^2 / d 2 k T ^ T (k T )^ T (k T -Q)^(Q,k T )FP(Q,k T ), 



where 



^(Q,k T ) 



= 1,2; M = 3,4,5) 



'exp(-^Q 2 ) 



fc = f, 



[exp { - [Q 2 " |Q • (2k T ) + |(2k T ) 2 ] }, k =/= I . 

For a proton-proton scattering, the spatial matrix element is 

J d 2 b e- lC > h J dq, d^dT7„ d£ n ¥£ f fo, £ x ) tfjy (??„, £ n )^T(rifeT)M^T) ^i.ifai, £i)*n,i(*7n> £n) 

= -^-a / d 2 k T VT(kT)VT(kT-QT)^(Q,kT)^(Q,k T ), = 1,2,3; fc,Z = 4,5,6) . 



(27T) 



(36) 



(37) 



(38) 



(39) 



V. TOTAL HADRON-HADRON CROSS SECTIONS 

After we have evaluated the color matrix element, the spin matrix element, and the spatial matrix element, we can 
now determined the total hadron-hadron cross section as 



where 



Ctot 



ijkl *ijkl ■) 



i,j,k,l 



(40) 



Cf,ijki = Cijki (Ci I A"A^ | Ci )( C u | A^Af | Ch ) 

,6=1 



(41) 



Hjkl 



(2tt) 
1 



l — J d 2 k T ^ T (k T )l/ T (k T )4 (Q-0,k T )^(Q = 0,k T ) 



d 2 k T < 



6nbn 



Ana* 



k 2 +^ 2 (k^+ZX 2 ) 2 



n 2 



(42) 



(21 

The superscript I (or II) is either M (meson) or P (proton), and ( s lfc j ; ) is the second-order spin-spin matrix element 



averaged over all polarization states of hadrons I and II, 

( S fk 3 l ) = JT (Xi Xn \( s i- s k){Sj ■ Si) | X : Xu ' 



(43) 



The polarization average of the first-order spin-spin matrix elements ( ) is zero. 

If we assume an unscreened linear-confining interaction for hadron-hadron collision (which corresponds to taking 
/i = in the 6irb / (kT + fj, 2 ) 2 term in Eq. (@2J0; the total cross section will be singular, since the integral over kp in Eq. 
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(1421) diverges. From a physical viewpoint, we expect the interaction between constituents in different hadrons to occur 
at relatively large separations. As a consequence, the interaction will be subject to screening, due to the production of 
virtual qq pairs. The linear confining interaction for bound states should thus become a screened confining potential 
when applied to the interactions of quarks in a diffractive hadron-hadron scattering process. Similarly we expect the 
color-Coulomb interaction to be screened. The spatial dependence of the interaction in Eq. (0 between quarks and 
antiquarks of different hadrons should be different from those inside the same hadron. 

A simple way to incorporate screening is to replace k 2 in the Fourier transforms of the color-Coulomb and the linear 
interactions by k 2 + p 2 , which leads respectively to the Yukawa and exponential potentials, as given in Eqs. © and 
(|7|) and the interaction in Eq. (|42|) . Previously, in studying the screened potential and its temperature dependence 
inferred from lattice gauge theory, the effective string tension bo = 0.35 GeV 2 was found to give good descriptions of 
meson bound state masses j^- We shall use this effective string tension in the present work. 

The rms radii of hadrons are model-dependent as the interactions are sometimes implicitly included in the determi- 
nation of the effective rms radius. We shall use the meson bound state wave functions obtained in the non-relativistic 
quark model of Ref. 0], where the rms radius is 0.256 fm for a pion and 0.261 fm for a kaon. (See Table IV of 
Ref. 01 1 m which the quantity \J (r 2 ) in the fourth column is the rms q-q separation, which is twice the rms meson 
radius.) In our present model in which a hadron-hadron reaction takes place via a finite-range interaction between 
constituents in one hadron with constituents of the other hadron, these hadron constituent radii can be considered 
"core" radii. These pion and kaon rms "core" radii are smaller than the corresponding values of 0.61 fm and 0.54 fm 
found in the geometrical model of Chou and Yang UllE} as weu as * ne vames 01 0-56 fm and 0.53 fm determined from 
electromagnetic measurements . These differences arise from the fact that Chou and Yang describe hadrons 

as a geometrical droplet with an essentially local zero-range interaction between a constituent of one hadron with the 
constituent of the other hadron. In the electromagnetic measurement, the photon also fluctuates into a p meson, and 
the larger electromagnetic rms radii include the additional effects of this fluctuation. 

It is reasonable to assume that the hadron radius niadron(EM), obtained in an electromagnetic measurement, is the 
sum of the hadron radius used here (the core radius) and an effective rho meson radius (which is of order l/m p ), 

^hadron (EM) = (hadron core radius) + (effective p meson radius). (44) 

Based on this assumption, the difference of of the proton and pion rms electromagnetic radii is equal to the difference 
of their corresponding rms (core) radii, 

^^(EM)) - y/(rl(EM)) = ^ ~ ^fW) ■ (45) 

The proton and pion rms radius obtained in electromagnetic measurements are 0.81 fm |50j and 0.56 fm [4^ | respec- 
tively. Using the meson bound states wave functions in Ref. 18] for which \J (r 2 ) = 0.256 fm, the above relation 

gives a proton rms radius of ^/(?' 2 ) = 0.51 fm which we shall adopt in our present analysis. This value of the proton 

rms radius is close to the value of «/ (r 2 ) = 0.48 fm found by Copley, Karl, and Obryk .51], and Koniuk and Isgur 

|52| in their quark-model description of baryon electromagnetic transition amplitudes. We shall use the quark masses 
and the hyperfine smearing parameter d that are employed in the meson bound state wave function calculations in 
Ref.0: m u = m d = 0.334 GeV, m s = 0.575 GeV, m c = 1.776 GeV, m b = 5.102 GeV, and d = 0.897 GeV. 

We search for the screening mass p and the strong coupling constant parameter a s by fitting the experimental cross 
sections of cr tot (p + p) = 39.0 mb, otot^+p) = 24.0 mb, and a to t{K + p) = 20.5 mb at ^/s — 20 GeV (53[. We obtain 
p = 0.425 GeV, and a s — 0.495 by using the Levenberg-Marquardt method p4j. 

The screening mass p = 0.425 GeV falls in between and m p studied by Dolejsi and Huffier j^]. It is interesting 
to note from Fig. 3 of their work that a screening mass between m v and m p gives the best correlation between the 
total cross section and the effective elastic scattering slope parameter b e ff. It is therefore important to extend the 
present calculations to the analysis of elastic scattering in future work. 

The screening length 1/p = 0.464 fm is the same physical quantity as the correlation length a studied by Landshoff 
and Nachtmann [31|. Using a single gluon propagator with a screening mass p, they noted that 1/p must be much 
smaller than the radius of a hadron, for the additive quark model to be valid. Our numerical analysis, with the full 
non-perturbative interaction of Eq. ©, leads to a range 1/p that is not small compared to the radius of a hadron; 
we should not expect our results to coincide completely with those of the additive quark model in all aspects. 

The value of a s = 0.495 is slightly smaller but close to the conventional strong interaction coupling constant of 
a s = 0.6 used in non-relativistic quark models for light quarks |21| . 

We shall first focus our attention on the total hadron-hadron cross sections at y^s = 20 GeV in Fig. 2 and then 
discuss the energy dependence in Eqs. (|51|) - (|5^|l and Fig. 3. The experimental pp, np, and Kp total cross sections 

2 



at yfs — 20 GeV can be roughly represented by a to t ~ 5.86(^/ (r 2 ) + Vvn)) 2 - While the theoretical total cross 
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sections obtained in the present TGMP calculations tend to increase with the square of the sum of the rms radii, 
there are however variations that arise from the spin-spin interaction. The latter quantity depends on the masses of 
the interacting quarks and the strength of the spin-spin interaction involving strange quarks is smaller than that of 
light quarks; consequently, ertot(?r + p) > &tot{K + p), as shown in Fig. 2. 

Our total cross section for p + p at y/s — 20 GeV is 27.0 mb, which is slightly larger than the total ir + p cross 
section of 24.5 mb. In the TGMP, the spin-spin interaction affects the p meson cross sections in two ways. First, 
as the quark-antiquark potential is weaker for the p meson than for the pion due to the spin-spin interaction, the p 
rms radius is 0.385 fm, which is substantially larger than the pion rms radius of 0.256 fm. Secondly, the spin-spin 
interaction also affects the interaction strength between a quark or antiquark in the p meson and a constituent in the 
other hadron. Hence, the TGMP predicts <r to t{p + p) slightly greater than <7tot(7r + p)- This result is slightly greater 
than the prediction of the additive quark model which gives 

o-tot(/0 + p) (additive quark model) = ^[crtot{^ + + P) + ^tot{^~ + p)] ~ 24.2 mb. (46) 

It is interesting to note that our TGMP total cross sections for ir+n, and p+ir are respectively 27.4 and 20.8 mb, 
which are significantly larger than the predictions of of about 2<7tot(?r + N)/3 ~ 16 mb from the additive quark model. 
The large magnitude of the total cross sections in the present calculation arises from the spin-spin interaction, as the 
cross sections become much smaller if the spin-spin interaction is turned off. 

Experimental information of the p+p cross section can be inferred from p- meson photoproduction cross sections. 
The extracted total cross section depends on the coupling constant /p/47r and the value of ry p , which is the ratio of 
the real and imaginary parts of the scattering amplitude |H3 • It also depends on the generalizations of the Vector 
Dominance Model in which one may include higher resonances and off-diagonal matrix elements |56J, |57j . 

Using the simple VDM with rj p assumed to be zero and /p/47r set to be 2.20, the ZEUS Collaboration found 

<Jtot{po+p) = 28.0 ± 1.2 mb (47) 

for y/s = 70 GeV t 58] which is close to the value of atot{p + p) = 30.5 mb obtained in the present calculation for 
y/s = 70 GeV (see Fig. 3). 

At low energies, the experimental data cannot be compared directly with the present theoretical prediction. One 
can nevertheless get some idea on the magnitude of a% Q t{p+p). Using photoproduction data from a number of different 
nuclei, McClellan et al. [5^ found earlier 

^tot(po+p) =38.0 ±3 mb (48) 

and /p/47r = 4.40 ± 0.60 for photons at 6 GeV {y/s = 3.48 GeV). However, from their subsequent extended measure- 
ments and refined analysis they obtain 

^tot(po+p) =25.9±1 mb (49) 

and /2/4tt = 2.52 ± 0.08 with rj = -0.24 60] for photons at 4.4-8.8 GeV {y/s = 3.02-4.17 GeVV 

From an extensive optical model analysis taking the ratio rj p to be -0.2, Alvensleben et al. |6l| found 

<rtot(pQ+p) =26.7±2mb (50) 

and /2/4tt = 4.40 ± 0.60 for photons at 4.6-7.2 GeV {y/s = 3.09-3.87 GeV). Other measurements give 

<Jtot(po+p) = 22.5 ± 2.7 mb (51) 

with /2/4tt = 2.52 ± 0.08 with rj = -0.24 [H for photons at 2.75-4.35 GeV {y/s = 2.44 - 3.01 GeV). 

An independent determination of the p — N total cross section, using the 7 + d — > p + d differential cross sections, 
was carried out by Anderson et al. who obtained 

o-tot(po +p) = 27.6 ± 0.6 mb (52) 

and /2/4tt = 2.80 ± 0.12 at a photon energy of 18 GeV {y/s =5.98 GeV) [63|. 

In our non-relativistic potential model, the p and lu have the same spin and spatial wave functions. Consequently, 
within the two-gluon model of the Pomeron, cr% t{p + p) = 0tot(k> + p). The value of CTtot( w + p) extracted from 
experimental photoproduction data depends again on the coupling constant f^/^ir and 77^. The extracted values of 
a tot {oj +p) range from 33.5 ± 5.5 mb jgj at 6.8 GeV {y/s = 3.694 GeV) to 25.4±2.7 mb at 8.3 GeV H3 {y/s =4.06 
GeV). 
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Our examination of the experimental p + p and u + p cross sections at low energies indicates that they appear 
to range from 23 mb to 38 mb and the uncertainties are large. The cross section depends on energy, and the only 
data point of er tot (p + p) at \fs = 70 GeV that can be direct compared with the TGMP model indicates approximate 
agreement between the theoretical result with experiment. However, much more work needs to be done to confront 
the predictions of the TGMP and the additive quark model with the results of photoproduction experiments. 

For the cf> + p reaction the present results give a theoretical total cross section of 17.1 mb at yfs = 20 GeV, which 
agrees with the naive additive quark model of 

<Jtot(<t> + p) = <?tot(K + + p) + a to t(K- +p)-Otot(7r _ +p)~17.0 mb. (53) 

The experimental photoproduction data give a cr tot (</> + p) which ranges from 9.2 mb to 17.6 mb for photons at 6.4-9.3 
GeV (y/s = 3.59 — 4.28 GeV), depending on the values of /J/47T and rj^ assumed. 

The TGMP gives a theoretical total cross section of (Ttot(J/ip+p) =7.6 mb at yfs = 20 GeV. From the absorption of 
J/ip in pA collisions, an effective Jjjp+p dissociation cross section of about 4-6 mb is used to described the dissociation 
process IS IS IS HIS El 113 El 113 • For the above analysis of J/ip production in Pb-Pb collisions at y/s = 17.3 GeV, 
the energy at which the J /if) or its precursor collides with a nucleon is probably about y/s ~8 GeV, corresponding to 
the collision of a produced J/ip nearly at rest with an incident nucleon in the nucleon-nucleon center-of-mass system. 
Thus, a theoretical result of Otot{J/ip + p) =7.6 mb in which the predominant component may be the dissociation 
cross section, is approximately consistent with the experimental dissociation cross section of about 4-6 mb. It remains 
necessary to determined more quantitatively how the J/ip + p cross section depends on energy, how much of the total 
cross section of J/ip +p can be attributed to dissociation of J/ip into an open-charm pair, and whether the observed 
absorption of J/ip inside a nucleus may involve an admixture of charmonium states with a color-octet component. 

The total J/ip+p cross section can also be extracted from the dipole model of photoproduction 0- I n this 
model, light cone wave functions for the virtual photon and charmonium are assumed, and a universal dipole cross 
section is introduced to fit proton structure function at small x and a wide range of Q 2 . The dipole model gives 
(Ttot(J/ip +p) — 4.4 mb at y/s = 20 GeV, which is lower than the value of 7.6 mb obtained here in the TGMP model. 

Experimental data indicate that at high energies total hadron-hadron cross sections increase slowly with energy. 
This energy dependence is beyond the scope of the present TGMP model. We can include such a dependence 
phcnomenologically, by assuming that the strength parameters of the interaction depend on energy. Specifically, we 
assume the energy dependence given by Donnachie and Landshoff 1] , and parametrize the energy dependence of the 
interaction strengths as 

a,(V5) = a s {V7s = 20 GeV)/(VS), (54) 



&o(V5) = Mv^ = 20 GeV)/(VS), (55) 

where 

y O 0.0808 , v -0.4525 

[/(v^)] 2 = ±i£ , (56) 



(X s 0808 + y s -0.4525- 



s/s=20 GeV 



X = 21.70 mb, Y — 56.08 mb, and s is in GeV 2 . The energy dependence of our parameters has been so chosen that 
theoretical pp cross sections follow the energy dependence of Donnachie and Landshoff 's parametrization of the total 
pp cross section, which represents a very good fit to the experimental pp data. 

Using such an energy-dependence, we extend our results to energies near y/s = 20 GeV and beyond. The lines 
in Fig. 3 show our results for total cross sections as a function of y/s. The symbols of circles, open-triangles, solid- 
triangles, open-squares, and solid-squares are experimental data for pp, n-p, ir+p, K-p, and K+p & respectively. 
It can be seen that in regions where data are available, our results are consistent with experimental data. Numerical 
results for total hadron-hadron cross sections are given in Table I. 

Recently, from the elastic and proton-dissociative p° photoproduction data obtained by the ZEUS Collaboration 
at HERA 66], the non-resonant contributions were analyzed in terms of the fluctuation of a photon into a 7r + 7r~ 
pair and their interaction with a nucleon |67| . A model-dependent analysis of the non-resonant contribution by 
Ryskin et al. [67^ . based on the Drell-Soding approach [68L l69l|. yields an effective total 7rp cross section of a(np) = 
31±2(stat)±3(syst) for a pion-proton center-of-mass energy of about 50 GeV [70(. The data of the ZEUS Collaboration 
is shown as the solid diamond point in Fig. 3 and is consistent with the present results of Fig. 3. 



12 



VI. TOTAL CROSS SECTIONS FOR MIXED-COLOR CHARMONIUM ON A PION OR A PROTON 

A charmonium system that is produced in a parton-parton collision may be in a nontrivial color state, out of 
which different pure color states may be projected pjrj l38j| . This initial state wave packet may acquire a transverse 
momentum and will travel through the nuclear medium. It is therefore clearly of interest to study the propagation of 
a colored cc in a color-singlet hadronic medium. 

The color wave function of a color-octet meson is 

i 3 

\° 8C ) = ^ E (c=l,2,...,8). (57) 

V Q,/3=l 

For the scattering of a color-octet meson with a color-singlet meson, we have 

and for the scattering of a color-octet meson with a proton we have 

'32/9 S cd , i=j, k = l, 

± ( Cf | AfA| | Cf ) ( Q\ | MXl \Ch) = \ Z^ 9 ^ 1 = 1 (59) 



,6=1 



2/9 5 cd , i^j, fc^i. 



The color state of a mixed-color meson can be expressed as 

8 

| c ,mxed- color ) = 'ft | C 1 ) + 78c | C*° ) . (60) 



c=l 

For a mixed-color meson scattering on a hadron, we have 



V ^ ( (J m i xe d— color | ^<?^ | (jrnixed— color ^ ^ £fl | | 1 



a, 6=1 
8 



l7i| 2 ( | A" A* | C{ ) + \lsc\ 2 (Ct \ I Q 6 

a, 6=1 L 

x(C I \ | AgA? | C& > 

Iti 



8c \ 



8 

12 r 



^{^^^^^(ciiA^ici; 

a, 6=1 

■ E l^l 2 f E ( PP C I A " A ? I C i C ) ( °ii I A * A ? I C n ) ) • ( 61 ) 



a, 6=1 



From Eqs. i|40H — l|42(l and (|61|l , the total cross section of a mixed-color meson scattering on a hadron can be 
expressed as 

a iot = (l-Ps)al+ X +P s a*+ l , (62) 

where P$ — X^c=i I T8c | 2 gives the total amount of color-octet mixing. 

In Fig. 4(a) and 4(b) we show the total scattering cross sections of a mixed-color charmonium on a pion or on a 
proton, as a function of the size of the cc pair and the amount of color mixing at y/s — 20 GeV. The wave function 
of the charmonium state is taken to be a single Gaussian wave function characterized by a root-mean square radius. 
The cross sections are insensitive to the spin state of the cc system. The results in Fig. 4 are the total cross sections 
for the collision of a spin-triplet S = 1 cc. It can be seen that the total cross section of a color-octet charmonium 
(P 8 = 1.0) scattering on a hadron remains finite even as the size of the charmonium approaches zero, and they change 
only weakly with the charmonium size. This is quite different from the case of a color-singlet charmonium (Pg, = 0.0) 
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scattering from a color-singlet hadron. For the case of the scattering of two color-singlet mesons, if we turn off the 
spin-spin interaction, the integrand in Eq. @ is proportional to 

{V T (r 1T - r 3T ) + Vr(r 2 T - i"4t) - ^ T (ri T - r 4T ) - Vr(r 2 T - r 3T )] 2 , (63) 

and for the case of the scattering of a color-octet meson with a color-singlet meson, the integrand is instead proportional 
to 

[Vr(riT - r 3T ) + U T (r 2T - r 4T ) - Vr(ri T - r 4 T) - Vt(*2T - *3t)] 2 

9 

+ -[F T (r 1T - r 3T ) - y T (riT - r 4T )][Ur(r 2T - r 3T ) - V T (r 2 T - r 4T )] . (64) 

We can see that expression (|64|) does not approach zero when the size of the color-octet meson approaches zero 
(i"it — * T2t), unlike expression H63JI . There are cancellations in the case of color-singlet hadron-hadron scattering, 
which are not present in the scattering of a color-octet charmonium state on color-singlet hadrons, as first pointed 
out by Dolejsi and Huffier |3^ . 

The color-octet cross sections for Pg, = 1-0 are approximately the same as those of Dolejsi and Hiifher |3^|. even 
though the potentials used are quite different. For the collision of the color-singlet charmonium (Pg = 0.0) on tt and 
p, the results at yj (r^ E ) — 0.2 fm in Fig. 4 are close to the results for J/ip + n and J/ip +p in Fig. 2 as the rms radius 
of J/ip in our calculation is 0.202 fm. 

VII. SUMMARY AND CONCLUSIONS 

We have calculated total hadron-hadron cross sections for the collisions of tt, K, p, 4>, D, J/tp, ip' , T, T', and the 
proton using the Low-Nussinov model of two-gluon exchange, as developed by Gunion, Soper, Landshoff, Nachtmann, 
Dolejsi, Huffier, and Wong. A set of meson wave functions obtained from a non-relativistic quark potential model 
is used in our calculations. A phenomenological potential including screened color-Coulomb, screened confining, and 
spin-spin interactions is employed to describe the gluon exchange between the constituents of the interacting hadrons. 
The screening mass and the strong coupling constant of the potential are obtained by fitting the total cross sections of 
p+p, tt +p, and K+p collisions at y^i = 20 GcV. We extend our results to higher energies using the phenomenological 
energy dependence obtained by Donnachie and Landshoff 0. 

We find that the total hadron-hadron cross sections correlate with the size of the scattering hadron system. There 
are, however, important variations arising from the spin-spin interaction. This leads to ert t(7r + p) > &tot(K +p), a 
larger p rms radius, and consequently a Ctot (p + p) that is slightly greater than <7 t ot (tt + p) ■ The experimental total 
tt + p cross section at 50 GeV obtained by the ZEUS Collaboration agrees with the theoretical result of the TGMP 
and the experimental ZEUS data point of the total p + p cross section at y/s = 70 GeV is close to the theoretical 
prediction. 

While the TGMP and the additive quark model give about the same cross sections in some cases, there are significant 
deviations because the effective spatial range of the quark-quark interaction in the present analysis, is not small 
compared to hadron radii. The total ir+ir, and n+p cross sections obtained here differ significantly from those of the 
additive quark model. Much more work remains to be done to confront the predictions of the TGMP and the additive 
quark model with experiments. It will be of great interest to obtain quantitative experimental measurements of these 
cross sections, using for example the production of pions in photon-photon collisions in which each photon fluctuates 
into a pair of pions, as in the Drell [68( and S6ding_|6|| approach in photon reaction processes. 

As was previously noted by Dolejsi and Hiifncr J32j, we find that the theoretical total cross section for a color-octet 
meson is quite large. The large cross section arises from the fact that in the interaction of a color-octet meson, 
the color interactions of both constituents in the meson interfere constructively, while in a color-singlet meson, the 
interactions interfere destructively. Therefore, a color-octet meson traveling through a color-singlet nuclear medium 
is expected to suffer frequent collisions and will lose a substantial amount of energy. These results for the color- 
octet cross sections may be used to place a constraint on the color-octet production amplitude of a heavy-quark pair 
produced in a nucleon-nucleon collision by studying its propagation in the nuclear medium. 

It will be of interest in future work to study the elastic differential cross section based on the TGMP, which can 
be compared directly with experiment. It would be useful to develop experimental and theoretical techniques to 
infer cross sections of unstable mesons on various targets, which would be valuable to discriminate between various 
theoretical models. 
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TABLE I: Numerical results for total hadron-hadron cross sections. 
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FIG. 1: A two-gluon exchange diagram contribution to the meson-baryon elastic scattering amplitude. 
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FIG. 2: Total hadron-hadron cross sections as a function of the square of the sum of the rms radii of the colliding hadrons. 
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FIG. 3: Total hadron-hadron cross sections as a function of energy. The symbols of circles, open-triangles, solid-triangles, 
open-squares, and solid-squares represent experimental data for pp, ir~p, vr + p, K~p, and K + p |53|l . The solid diamond point 
is from the n+p data of the ZEUS Collaboration |6fij . 




FIG. 4: (a) The total cross sections for mixed-color charmonium scattering on a pion as a function of the charmonium rms 
radius at ^/s — 20 GeV. Ps specifies the probability of the charmonium color-octet state, (b) The total cross sections of 
mixed-color charmonium scattering on a proton as a function of the charmonium rms radius at ^fs = 20 GeV. 



